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In Brans-Dicke theory a non-linear self interaction of a scalar field (f> allows a possibility of real- 
izing the late-time cosmic acceleration, while recovering the General Relativistic behavior at early 
cosmological epochs. We extend this to more general modified gravitational theories in which a de 
Sitter solution for dark energy exists without using a field potential. We derive a condition for the 
stability of the de Sitter point and study the background cosmological dynamics of such theories. 
We also restrict the allowed region of model parameters from the demand for the avoidance of ghosts 
and instabilities. A peculiar evolution of the field propagation speed allows us to distinguish those 
<«_^ ' theories from the ACDM model. 

(N' 

^ ■ I. INTRODUCTION 

I— J " 

^T ' Scalar-tensor gravitational theories have been widely studied as an alternative to General Relativity (GR). In these 

theories a scalar- field degree of freedom (f) is coupled to the Ricci scalar R through a coupling of the form F{(f))R 
[l| . For example, such a coupling arises in low energy effective string theory as a result of the dilaton coupling with 
^^ . gravitons [2[. The discovery of dark energy in 1998 Q has also stimulated the study for the modifications of gravity 
r\ ' on large distances (see Refs. [j] for recent reviews). 

• , The well-known example of scalar-tensor theories is Brans-Dicke (BD) theory Q in which F(0) is proportional 

1^ ■ to (p with a non-canonical kinetic term (— a;BD/'/')(V(/))^, where ojbd is the so-called BD parameter. In original BD 

II [ theory without a field potential, the BD parameter is constrained to be wbd > 40000 from local gravity tests in the 

O ■ solar system [6|. This comes from the fact that the coupling between the massless field (f) and non- relativistic matter 

|H I needs to be suppressed to avoid the propagation of the fifth force. Under this bound, the deviation from GR (which 

^ ■ corresponds to the limit wed — >■ oo) is too small to be detected in current cosmological observations. 

I i ' There are two ways to recover the General Relativistic behavior in high density regimes relevant to solar system 

experiments. One is the so-called chameleon mechanism Q in which the presence of the field potential V{4>) allows 

CN ' a possibility of having a density-dependent effective mass of the field. Provided that the effective mass is sufficiently 

large in the regions of high density, a spherically symmetric body can have a thin-shell around its surface so that 

the effective coupling between (j) and matter is suppressed outside the body. This mechanism works not only for BD 

theory [3,H but also for metric f{R) gravity Q, because the latter can be regarded as a special case of BD theory with 

^^ ■ wbd = [i3, [ill- For viable dark energy models based on BD theory and metric f{R) gravity the early cosmological 

evolution mimics that of GR, but the deviation from GR becomes important at late times (i.e. low density regimes) 

^^ ' Another way for the recovery of GR in the regions of high density is to introduce non-linear self interactions of a 

scalar field, e.g., (,{4>)0(j){d^ (j)d^(j)) , where C is a function in terms of (p. There have been attempts to restrict the form 
of the Lagrangian by imposing the "Galilean" symmetry dfj,(f> — )■ dfj,(f> -\- h^ [1Jt13- The self interaction of the form 
D(t){d'^(l)d^(j)) respects the Galilean symmetry in the Minkowski background. In the Dvali-Gabadadze-Porrati (DGP) 

S^ i braneworld model [l^ the field self interaction 0(f>{d''''(f>dfj,(f>) has been employed to recover the General Relativistic 
5—1 ' behavior for the length scale smaller than the so-called Vainshtein radius [20, UM, (see also Refs. [22|). While this 
"self screening" mechanism (called the Vainshtein mechanism [23] ) can be at work for consistency with solar system 
experiments, the DGP model is unfortunately plagued by a ghost problem [2l|, [2J] as well as incompatibility with a 
number of observational constraints [2g ■ 

In BD theory with the field self interaction term £^{(f>)0(f>{d'^ (pd^cj)) , there exist de Sitter (dS) solutions for ^(0) ex (f)^'^ 
and wbd < —4/3 ^^^ . Moreover the General Relativistic behavior is recovered in early cosmological epochs during 
which the field is nearly frozen. The solutions finally approach the dS attractor at which (p/cf) is constant. As long as 
(j)/(f> is positive, one can avoid the appearance of ghosts and instabilities associated with the field propagation speed. 
Moreover this model gives rise to several interesting observational signatures such as the modified growth of matter 



> 

00 

00 



o 



a 



^ For constant ^{<p) there exist no consistent dS solutions, but the presence of other field self interaction terms that respect Galilean 
symmetry in the Minkowski background allows a possibility for giving rise to dS solutions [18J . 



perturbations and anti-correlations in the cross-correlation of large scale structure and the integrated Sachs- Wolfe 
effect in cosmic microwave background anisotropies [2^, [23| • 

In this paper we consider the general action ([1]) below without the field potential and derive the functional forms 
of F((f)), B{(j))^ and ^(0) from the requirement of obtaining dS solutions with (j)/(j) = constant. The three functions are 
restricted to be of the power-law forms in terms of 0, which include BD theory with ^((/)) ex ^~^ as a special case. We 
discuss the cosmological viability of such theories by analyzing the stability of fixed points for the background field 
equations. We also study conditions for the avoidance of ghosts and instabilities to find the region of viable model 
parameters. 



II. GENERALIZED BRANS-DICKE THEORIES 



+ / d4x£M(5M-,*A/), (1) 



We start with the action 

S = / d^x. 



\f{cj,)R + Bi<j>)X + a^)D^id''^d^^) 

where g is a, determinant of the metric g^^^,, (/) is a scalar field with a kinetic term X = —g^^'^dp,(l)d^(j)/2, and F{(f)), 
-S(^), C(^) are functions of 0. Cm is a matter Lagrangian that depends on the metric g^jy and matter fields ^a/- We 
would like to construct models in which the late-time cosmic acceleration can be realized without the field potential. 
We shall restrict the functional forms of i^(0), 5(0), and ^{(j)) from the requirement of having dS solutions. In the 
presence of the nonlinear field self-interaction term it is possible to recover the General Relativistic behavior for the 
cosmological evolution at early epochs. 

We study the cosmological dynamics for a spatially flat Friedmann-Lemaitre- Robertson- Walker (FLRW) space-time 
with the line element 

ds^ = -dt^ + a{tfS^■i dx'dx^ , (2) 

where a(t) is the scale factor with cosmic time t. For the matter Lagrangian Cm we consider perfect fluids with 
non-rclativistic matter (energy density pm) and radiation (energy density pr and pressure Pr = PrI'i)- Varying the 
action ([Ij with respect to g^y and </> for the metric ([2]), we obtain the following equations 

3Fif2 = B0V2 - 3ffF^0 + {mi - C,00)0^ + Pra + Pr , (3) 

-2FH = {F^^ - 2e</)2)0 + [B0 + ^,000 - HF^^ + <j)^i6H( - 2i^<P)]<p + p,„ + 4p,,/3 , (4) 

{B - 402^,^ + 12HC^)4> + [3HB + B^^^/2 ~ i^^<j>^ + 6^i3H^ + H)^]^ ~ 3F42H^ + H) = , (5) 

where a dot represents a derivative with respect to t, and H = a/a is the Hubble parameter. Note that we have used 
the notations F^ = dF/d(j) and F^^^ = d'^F/dcfP' (the same for other variables). 
From Eqs. (O and dj) it follows that 
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2 (P^-<t' 



^ + f^i/2^.3 1 _ :^x + -^ + -^ , (6) 




-24=r^-^i.vyA+.A+.,.U[^x-f^>^.-^-K^i7vri-^. 



where 



H^ \ F F ' J \H H'^ ' j [ F F F F \ ^ 

Pm 4:pr 

FH'^ SFH^ 



H(j) 



(7) 



(8) 



Now we look for dS solutions at which H and x are constants. One can generalize the analysis to the case in which x 
is a function of 0, see the Appendix A. Here we focus on theories giving constant x, to recover the dS solution found 
in Ref. [2y| as a specific case. If F{4>) and ^((/)) are power-law functions in terms of 0, the quantities such as cj)F^^/F, 
F^^^(f>^ /F, and 4i£,.,f,/£, remain constants. Provided that B/F oc 0~^ and £^/F ex 0"'^, one can solve Eqs. ([6|) and ^ 
for x and H at the dS point. From the above demands we adopt the following functions 

F(0) = M2i(0/Mpi)3-" , B{<j>) = c.(0/Afpi)i-" , m = (A/M')(0/Afpi)-" , (9) 



where Mp\ ~ 10^^ GcV is the reduced Planck mass, ^ (> 0) is a constant having a dimension of mass, and uj and A are 
dimensionless constants. For the consistency of theories, the coupUng A needs to be positive. In the Appendix B we 
shall show why the theories with A < are not allowed. The Brans-Dicke (BD) theory [5| corresponds to n = 2 with 
the BD parameter uj. Since F{(f)) is constant for ti = 3, the theory with ti = 3 corresponds to k-essence [29| minimally 
coupled gravity. Note that the choice of functional forms different from ([9|) may allow the possibility to realize the 
cosmic acceleration today, but we shall focus on the simplest case in which dS solutions are the late-time attractor. 
From Eqs. ^ and ([7]) we obtain the following algebraic equations at the dS point: 

n(n-3)2x3g + (n-3)(n-12)x2g-6(n-5)xds + 18 ,^„. 

a;ds(a;ds + 3) 

^ ^ M^ [{n - 3)a;ds - 2] [(n - 3)a;ds - 3] 

A/piH^g 2.T3s(xds+3) ' ^ > 

where Xds and H^s are the values of x and H at the dS point, respectively. We fix the mass scale /x to be 

Pi = {MpiHis)'^^ ^ 10-^Hlpu (12) 

where we have used H^s — 10~^°Afpi. In this case Eq. ([TT]) is simplified as 

[(n - 3) xds -2][(n-3 )xds - 3] 
2a;3s(a;ds + 3) 



> ^ LV"- - oy-i-dB - ^JLV"- - ^J-^dii - ^l ,^„. 



which gives Xds ~ C'(l) for A and n of the order of unity. For given uj and n, the quantity Xds is determined by 
solving Eq. pU)) . Then the dimensionless constant A is known from Eq. p^ . The demand for realizing the cosmic 
acceleration today fixes the value of A together with the mass scale fi given in Eq. ((T2|) . 
For BD theory (n = 2), Eq. pU)) can be analytically solved as [2a| 

^2±V-8-6uj 

a;ds = TT- . -3 ■ (14 

2 + UJ 

As we will see later, the condition for the avoidance of ghosts demands that x > during the cosmological evolution 
from the radiation era to the dS epoch. When n = 2. we obtain x^g > only for the solution a;ds ~ (—2 — 
V-8-6w)/(2 + w) with UJ < -2. 

In order to recover the General Relativistic behavior in the early cosmological epoch we require that the field initial 
value 0i is close to Mp\. As we will see later, the quantity x is much smaller than 1 in the early cosmological epoch. 
Hence the field is nearly frozen during the radiation and matter eras. The field starts to evolve at the late cosmological 
epoch in which x grows to the order of unity. 

For the functions (|9]), let us derive autonomous equations of the dynamical system. In addition to the quantity x 
defined in Eq. (j8]) we also introduce the following variables 

,^A.2^, ^r^^. (15) 

Since we are considering the theories with A > 0, the variable y is positive definite. Equation (|6]) gives the constraint 
equation 

"'" ^ sl^ = 1 - f^r - f^DE , (16) 



where 

n - ^ 
"DE = -zx 

6 
The variables x, y, and 51^ obey the following differential equations 



(3 - n)x + 2xy (l + ^x^ . (17) 



d2^ 9x ^^ ^ n.^fr, , ^n,\ , 6(3 - n)a; 2 ^H' 



-x~ -—(3 -n-2y)[n„, + -17, + ^ , ' ~ x' - x— , (18) 

di Zujdi \ 3 / ujdi rl 



y' = -2y 



d2 9 ,^ ^ , / 4^ \ 6(3 -n) 

(3-?i-2y) a« + ^a- - ^ , ' +x 



di 2ujdi \ 3 / ujdi 



(19) 



n',^-2nr(2+'^x + ^] , (20) 



where a prime represents a derivative with respect to TV = In a, and 



d2^ 


12 

= — y + X 

UJ 

18 

= — y + a;< 


l + -y+^(3-n-2j;)2 

id ZLO 

1 — n nin + 1) 


5 

3(3 -n- 2y) 


2 a; 



X (3 — n)(2 — n) 3 — n 
2 21^; ^ 2^ 



-2; 1 



(21) 
(22) 



The Hubble parameter satisfies the equation 
H' 1 



H 



:{i-n~2y) 



d2 9x ,„ „ , (^ 4^ \ 6(3 -n) 
X + ::: — —\o — n — 2y) [ll„i + —\lr I ; x 



di 2ujdi 



Cddl 



1 3 

--X [u!x + (3 — n)(2 — n)x — (3 — n) + j/(6 + 2nx)] — -Sl,„ — 2$!,. 



We also define the effective equation of state 



Wcff = -1 - 



2H' 



(23) 



(24) 



The relation ([23]) should be substituted into Eqs. ([T5|) . (pn|) . and ([M|) to solve the cosmological dynamics numerically. 
At the dS point (H = constant) we have that Woff = ^1- 

In addition to the dS point derived above, there is a fixed point for the system ([T5| - (pn)) that corresponds to the 
matter-dominated epoch: 



P^: ix,y,nr)= (0,^—^,0 



(25) 



Note that this can be also derived by equating the term 6^(3i?^ + H)<i)^ in Eq. ^ with the term 3F^^{2H^ + H). 
For A and n of the order of unity we have that Xds = ^'(1). From the definition of y given in Eq. (jisp we have 
x^ ~ y{H^a,/ H^) <C y during the matter and the radiation eras. Since y is positive definite, we require the condition 



n< 3. 



(26) 



Considering linear perturbations 5x^ Sy, 6 fir about the fixed point (j25p . the eigenvalues of the Jacobian matrix of 
perturbations are given by 3/2, —3/2, —1/2. Hence the matter point is a saddle followed by the dS solution. 

There are no fixed points that correspond to the radiation-dominated epoch. However one can analytically estimate 
the evolution of the variable y under the condition that x is negligibly small relative to y. Using the approximation 
that di ~ 12y/uj, d2 — l&y/uj, and 17^ ~ 1 — ^m, we obtain the following equation during radiation and deep matter 
eras: 



y' - -y + t(3 - n - 2y)nm ■ 
Provided that 2y <^ 3 — n, the solution to Eq. ()27p during the radiation domination is given by 



(27) 



y- -{i - n)n,n + ce' 



-N 



(28) 



where c is a constant. To derive the solution (|28p we have used the relation ft'^ ~ 51„i ex a. Since the term ce~^ 
decays with time, one has y ~ (3 — n)i7,„/8 ex a for 2y <^ 3 — n and n 7^ 3. When the term 2y grows to the same order 
as 3 — n during the matter dominance (£7^ — 1) one can obtain the matter point (|25p by setting y' = in Eq. (j27p . 
Note that if n = 3 the variable y decreases as y ex e~^ during the radiation domination. 

The above discussion shows that the radiation era, characterized by (x, y,^r) — (Oj (3 — n)i7m/8, 1) and Weff = 1/3, 
is followed by the matter era with {x,y,ilr) — (0, (3 — n)/6,0) and uicff = 0. The solutions finally approach the dS 
point characterized by the conditions (|10p and ((TT|) . provided that it is stable. In the next section we derive the 
stability condition for the dS point by considering homogeneous perturbations about it. We also restrict the viable 
region of model parameters from the requirement to avoid ghosts and instabilities. 



III. STABILITY OF THE DE SITTER SOLUTION, GHOSTS, AND INSTABILITIES 

In order to discuss the stability of the dS solution as well as conditions for the avoidance of ghosts and instabilities, 
we consider scalar metric perturbations a, /3, ip, and 7 about the flat FLRW metric [30| 

ds^ = -(1 + 2a)dt^ - 2a{t)d,(3 dtdx' + a^{t) (% - 2^:6^^ + 2d,dj-f) dzMx^' . (29) 

We introduce the gauge-invariant curvature perturbation [3l| 

TT 

n = i' + ^s<j) , (30) 

which can be used to discuss the stability of cosmological solutions. 

A. Stability of the dS solution 

We expand the action ([T} at second order in the perturbation around the dS background. There is only one 
propagating scalar degree of freedom, which corresponds to the curvature perturbation TZ. In terms of the gauge 
invariant quantity TZ, the second-order perturbed action about the dS point is given by 



(55(2) ^ I ^td^xa'^Q, 






(31) 



where 



2Fr 
Q^ ^ ,^ ; . ^„ ^ ^.-.... ^ (32) 



cl- ' 






r(F^ - 2^<i?) 

-72HdsFF^^<j>^ - 3(24^ij2gi^2 _ p3.^2 ^ uh^^FFU + 12H%F^F^ 



(33) 



and 



12^^06 „ 4(3^^^^ ^ ^^^^)^4 ^ ^p2j2 ^ i2H^^FF^<i) + 12H%F^ . (34) 



The scale factor evolves as a ex cxp(iJdS t) at the dS point. If ^ = 0, then the propagation speed reduces to unity. 
For the power-law functions of F[(f)), B{(j)), and ^((/)) given in Eq. ^, Eqs. ([5^ and ([55)1 reduce to 



2k1-"(/)3-« [(„ _ 3)(3^2 _ iQ^ ^ I2)xlg + 12{n - 5)4g + 18(n - 8)a;ds - 108] 

(n - 2)2 [{n - 3)xds - 2] x^g 
{n - 2)[{n - 3)(n - 4)^3 ~ 8(n - 3)a;ds + 6]a;ds 
(n - 3)(3n2 - lOn + I2)xlg + 12{n - 5)xlg + 18(n - 8)a;ds - 108 



(35) 

yn — ^j- yyii — OjXdS — ^J ^dS 

2 _ (n - 2)[(n - 3)(n - A)x\^ ~ 8(n - 3)a:ds + 6]a:ds ^ocn 



From the action pip we obtain the equation for the curvature perturbation in Fourier space: 

^ '^-{a''Q,n) + cl%n^O, (37) 



a^Qs dt 



where A; is a comoving wavenumber. For homogeneous perturbations {k = 0) which have only the time-dependence, 
the solution to Eq. ([57)1 is given by 

n[t) = ci + C2 / -^ dt , (38) 

where ci and C2 are integration constants. Since x^s is constant at the dS point one has <f> oc exp(a;ds^dsO and hence 
Qs oc exp[(3 — n)xdsHdst] from Eq. (j35p . Then the homogeneous perturbation about the dS point evolves as 

n{t) = ci + C2 exp{[(n - 3)a;ds - 3]Hds t} , (39) 



where £2 is a constant. In order to avoid the growth of 7?,, one requires that 

{n ~ 3)xds < 3 , (40) 

which corresponds to the stabihty condition for the dS point. Note that the same condition follows by considering 
the stability of perturbations 5x, Sy, and 6flr for Eqs. (|18p - ((20|) . Since n < 3, Eq. ([40|) is automatically satisfied for 
s^ds > 0. In fact, the variable x needs to be positive during the cosmological evolution starting from the radiation era 
to avoid the appearance of ghosts. In the following we will show how this condition arises. 

B. Conditions for the avoidance of ghosts and instabilities 

Let us derive conditions for the avoidance of ghosts and instabilities by taking into account non-relativistic matter 
(equation of state Wm — +0) and radiation (equation of state Wr ~ 1/3). The velocity potentials «,„ and Vr are 
related with the energy-momentum tensors T^ and T? of non-relativistic matter and radiation, respectively, as 
Tj? = —PmVm.i and Tf ~ -^{Pr + Pr)vr.i- We introduce the following gauge-invariant combinations 

Scf) d<f> 

Vm=Vra ^ , K- = "r j- ■ (41) 



Perturbing the action ([T]) at second order, it follows that 



(55(2) = ! / dtd3.Ta3 



V^AV - VV^PW - V^BV - V^MV 



(42) 



where V = {TZ, Vm, Vr), and A, P, B, M are time-dependent 3x3 matrices. Here we do not write the explicit forms 
of the matrices because of their complexities. The readers may refer to the papers [33, l3j] for the details of such 
studies. 

The ghost conditions can be derived by the signs of the eigenvalues of A. Diagonalizing the 3x3 matrix A with 
components Aij, we find that the ghosts are absent under the conditions 

detA>0, A22A33 - Als > , A33>0. (43) 

In the uniform field gauge given by S(f> — 0, these conditions translate into 

dct A - P-^^ + ^-) ^-(^ + ^'■) ^-^'^'[^^.l + ^^^ + ^^^'^' - 4(2-^^^ + g^^^O-^^ + 24gf e<^] ^ p 
Wm Wr [F^^4> + 2FH~2^4y^Y 

A,,A,, A\, = Pma + -^r.)Pril+Wr) ^ ^ ^^^^ 

Wm Wr 

Pr(l + Wr) 

A33 = — > . 46 

Wr 

Since the conditions P5|) and (pS)) automatically hold, the no-ghost condition is simply given by det A > 0, i.e., 

2F02[3f2 + 2FB + 12^204 _ 4(2^^,0 + 3F^O0^ + 24i7Fe(/'] 



2FH - 2^(/)3) 



3^2 



> . (47) 



Eliminating the term B by using the background equation ([3]), the above definition of Qs is identical to (|32| at the 
dS point. For the power-law functions given in Ec^. (|9]) this condition reduces to 

2Fx[x{27-18n + 3n^ + 2uj + 12y^ + 2Q7iy~36y) + 2Ay] 
^'~ Kn-3 + 2y)-2p >°- ^^^> 

During the radiation and matter eras one has |a;| ^ y < 0(1) and hence 

Q, ~ 12i^a;y > -^ Fx > . (49) 



The no-ghost condition arising from vector and tensor perturbations corresponds to i^ > 0. We then find that Eq 
gives X = (j)/{H(f)) > 0. As long a.s cf> > initially the field and the function F{(j)) = Mpj((^/Mpi)'^~" grow for n < 3 



which means that F{(j)) remains positive. If the variable x crosses from positive to negative, this corresponds to the 
violation of the no-ghost condition because Qs becomes negative just after the crossing. Hence we require that .t > 
during the cosmological evolution starting from the radiation era. Under the condition (|26p together with F > and 
X > 0, we find that Qg in Eq. psp is positive at the dS point. 

Let us next discuss the three speeds of propagation. Two of them are trivial, c^„ = Wm , c^ = Wr, whereas the third 
one can be written as 

cl = ^2 (-48?'*'/''' + 64F?2^,^(/)i° ~ 384i7F?3^9 + [(72^2 + ASFF,^^)^ + (-32FF^C,0 + ^QF^i^^)^ 



Z2 



^16F^^%]4>'' + {-MHF^^t4> + 576HFF^^e)4>'' + [{672F^H^ + ((48uv„ - 80)p™ + (48u;, - 80)pr)F)e 
-i24FF^^F^^^ + 8F^B,^ + i8F^^)^]<j)'^ - 336HFF^^^<p'' + [{-624H^F^F,^ + ((88 - 24w,„>„, 
-(88 - 24wr)pr)FF^^)^ + 9F^^]^^ + [{-38iH^F^ + 128(p™ + Pr)HF^)£, + 72HFF%]4,'^ 
-\216H^F^F^^ - 24(p,„ + Pr)FF^^]^^ + [288H^F^F^^ - 96(p,„ + pr)HF^F^^]<j) + lAAH^F^ 

-96{p^ + Pr)H^F^ + 16{pr,^ + Pr)^ F^} , (50) 



where 



Z=(P^ [3F' + 2FB + 12e^0'* - 4(2i^^,^ 4- 3F^^C)r + 2AHFC(j)] . (51) 



We require that c^ > to avoid the instabilities for small-scale perturbations. 

Plugging Wjn = and Wr = 1/3 into Ec^. ([SH]) and using ilr and il,„ defined in Eqs. ((TS|) and ((T5)) . it follows that 

cl = {[-48/ + Wny^ + (-82477. + I92n'^ + 936 + 40uj)y^ + (-1728 + 68nuj - 140w - 76871^ + 9677.^ + 2016n)y 
- 72nuj + Au? - lOSn^ + 12^^^ + 729 - 972n + 108w + 48677^ + 9n^\x^ + [96y^ + (-576 -F 3207i)y2 
-f (32a; + 648 - 43277 + 72772)y]a; + 48(17^ -f h)y^ + 24(3 - 7i)(l - ^r)v] 
X [(12y2 + (2071 - 36)y + 377^ + 2w + 27 - 1877> + 2AyY'^ . (52) 

Let us estimate the evolution of c^ during the radiation and deep matter eras in which x is negligibly small relative 
to y. Using the relation fJ^ — 1 — f^m, we obtain 

1 (3-n_2y)a^^ 
" 2 24y ^ ' 

During the radiation domination the variable y evolves as y ~ (3 — 77,)0„i/8 for 77 ^ 3, see Eq. (|28l). Plugging this 
relation into Eq. (1531) . we find 

5 1 5 

c^ ~ r^m « - , (radiation era for 77 7^ 3) , (54) 

6 12 6 

which is valid for ri,„ ^ 1. During the matter dominance the evolution of y is given by y = (3 — ?i)/6. From Eq. (j53p 
we then obtain 

2 
c^ ~ - , (matter era for n 7^ 3) , (55) 

which agrees with the result derived for 77 = 2 [23. This propagation speed is independent of the values of 77. As 
we will see in the next section, c^ starts to deviate from 2/3 around the late epoch of the matter era. The above 
discussion shows that the instabilities associated with the negative propagation speed squared can be avoided during 
the radiation and matter eras. 

While the results (|54p and ([55|) have been derived for ri 7^ 3, the situation is different for n = 3. Plugging n = 3 
into Eq. (|53p , the propagation speed squared in the regime a; ^ y is given by 

c^ ~ - — — fim , (early cosmological epoch for n = 3) . (56) 

As we see in the next section, this estimation ceases to be valid during the matter era because x and y become the 
same order for n = 3. 
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Figure 1: The propagation speed squared at the dS point versus Xds for several difTerent values of n. One has < c^ < 1 for 
2 < n < 3. 



At the dS fixed point, c^ is given by Eq. (p6| . Under the conditions n < 3 and Xds > 0, one has c^ > for 

n> 2. 



(57) 



In Fig. [T] we plot c^ at the dS point as a function of Xds for a number of different values of n. One has c^ = for n = 2. 
In the regime Xds ^ 1 the propagation speed squared is approximately given by c^ ~ (n — 2)a;ds/18. If 2 < n < 3 the 
asymptotic value of c^ in the limit Xds — >■ oo is 



(n-2)(4-n) 



lim c„ — .-, , 

:jds^oo " 3n^- lOn+12 



(58) 



which has the maximum value c^ = 0.134 at n = 2.618. In Fig. [T] we find that, for n < 2.9, c^ monotonically 
increases with Xds- The local extrema of c^ appear for n close to 3. When n = 3 the propagation speed squared has 



a dependence cj ~ l/(4xds) for Xds 3> 1 and hence cj approaches in the hmit that Xds 
discussion we have < c^ < 0.134 for 2 < n < 3. 



From the above 



C. Viable parameter space 



Let us summarize the viable region of model parameters. From Eqs. (j26p and (j57[) the parameter n is restricted in 
the range 



2 < 71 < 3. 



(59) 



The condition 2 < n comes from the requirement c^ > at the dS point, whereas another condition n < 3 is needed 
for the existence of the matter fixed point. Provided F > 0, we require that 



x> 0, 



(60) 



to avoid the appearance of ghosts during the cosmological evolution from the radiation era to the dS epoch. Under 
([59)1 and (|5D|) the condition (^0)) for the stability of the dS point is automatically satisfied. 
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Figure 2: The allowed and excluded regions in the {n,uj) plane. The vertical axis is plotted in terms of logjQ(— oj) with oj < 0. 
The allowed region corresponds to 2 < n < 3 and uj < —n{n — 3)^. The black line shows uj = —n{n — 3)^, whereas the dotted 



line correspond to the border at which the propagation speed squared temporally reaches c^ 
cosmological evolution. 



1 during the course of the 



From Eq. (|10p one can restrict the parameter region of uj under the conditions (j59p and (j60p . In the limit that 
Xds ^ oo the asymptotic value of uj is lUoo ~ ^n{n — 3)^. Meanwhile, in the limit that Xds ~^ Oj we have uj ~ 
-CX3. The viable dS point exists for 



'6/^Is 



Lu < —n{n — 3) 



(61) 



vifhich does not have a lower bound. When n = 2 and n — i the allowed region corresponds to w < —2 and w < 0, 
respectively. In Fig. [5] we illustrate the viable parameter space in the (n, w) plane. The black line represents the 
border characterized by a; = —nin — 3)^. 

Although we have shown that c^ exists in the regime < c^ < 1 during the radiation, matter, and dS epochs under 
the conditions (|59|) and (j60p , it can happen that the speed of propagation temporally reaches the superluminal regime 
(Cj > 1) during the transition from the matter era to the dS epoch. We shall discuss about this in the next section. 



IV. COSMOLOGICAL EVOLUTION 



In this section we shall numerically integrate the dynamical equations (|18p - (j20p to confirm the analytic estimation 
given in previous sections. When n 7^ 3 the variable y evolves as y ~ (3 — n)r2,„/8 ex a during the radiation era, which 
is followed by the matter era with y ~ (3 — ri)/6 = constant. Finally the solutions approach the dS point with 



yds 



[{n - 3)xds - 2][(rt - 3)xds - 3] 
2a;ds(a;ds +3) 



(62) 



For n and lo satisfying the conditions ([5^]) and (pT|) . Xds (> 0) is known by solving Eq. ((TU)) . 

In Fig. [3] we plot the evolution of r^DE, ^m, and $7^ as well as x and y versus the redshift z for n = 2.5 and uj ~ —10. 
In this case Eqs. (fTO)) and ((62|) give Xds = 0.974 and yds = 1.120, which coincide with the numerical values at the 
dS point. The numerical evolution of the variable y shown in the right panel of Fig. [3] is consistent with the analytic 
estimation of y for each cosmological epoch. The evolution of the variable x during radiation and matter eras can be 
understood as follows. Using the relations x -^ y < C(l), O,. ~ 1 — Sl,„, and H' / H ~ — 3rj„i/2 — 251^ in these epochs, 
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Figure 3: Evolution of Ode, f^m, fir, and Weff versus the redshift z = 1/a — 1 for n — 2.5 and a; = — 10 (left panel). The 
initial conditions are chosen to be a; = 10~^* and y — {3 — n)Qrn/8 with Qm = 1-28 x 10~^. We identify the present epoch 
(z = 0) as f^DE = 0.72, Qrn = 0.28, and £7^ = 8 x 10~^. The right panel shows the evolution of the variables x, y and 17,. with 
a logarithmic scale for the same model parameters and initial conditions as those in the left panel. 



Eq. p^ reduces to 



3 (3 ~ n)fl„i ^ 3 
2 + Sy 4"^ 



(63) 



Since y ~ (3 — 7i)f2„i/8 for n 7^ 3 during the radiation dominance, we obtain the solution x oc e^^" = a^". During 
the matter dominance characterized by y ~ (3 — n)/6 and flm — 1, the solution is given by a; ex e^^/^ = a^^^. Hence 
the variable x evolves faster than y, as we see in the right panel of Fig. [3l 

The left panel of Fig. [3] shows that the successful sequence of radiation {wcb — 1/3, fir — 1), matter (woff — 0, 
Clm — 1); and dS eras (wcn — —1, SIde — 1) arc in fact realized. Even if wc start integrating from the high- redshift 
regime the solutions arc not prone to numerical instabilities, unlike modified gravity models having a large mass in 
the regions of high density (such as f(R) dark energy models [l3,l3J, |35|). This property is associated with the fact 
that, instead of the field potential, the scalar-field self interaction is used to recover the General Relativistic behavior 
at early times. 

We have run numerical simulations for other model parameters satisfying the conditions ()59p and (|6ip and confirmed 
that the successful cosmological evolution can be realized at the background level. We have also followed the evolution 
of the quantity F = Mp[(0/Mpi)'^~" by numerically solving the equation 



F' ^{3~ n)Fx . 



(64) 



Provided F > Q and x > initially, the numerical simulations show that F continues to increase for 2 < 71 < 3 
without changing the sign of x. In the simulation of Fig. [31 for example, the present value of F is about 1.35 times 
as large as its initial value at the radiation era. This growth of F is different from f{R) dark energy models in which 
the quantity F = df/dR decreases with time [IJ, ll3l- Numerically we have also checked that the no-ghost condition 
Qs > is satisfied from the radiation era to the dS epoch. When n = 2. Eq. ((35l) shows that Qs ^ 00 at the dS 
point. Hence the model proposed in Ref. [26| may suffer from the lack of proper quantizations at the dS point. 

In Fig. m we plot the evolution of c^ for several different cases. When n 7^ 3 the numerical simulations confirm the 
values C3 ~ 5/6 and c^ ~ 2/3 during radiation and matter eras, respectively. Depending on the values of n and w, 
the evolution of c^ from the matter era to the dS epoch is different. When n = 2.5, for example, the propagation 
speed remains subluminal (c^ < 1) for w > —85 (as in the case (b) in Fig. U), whereas it temporally reaches the 
supcrluminal regime (c^ > 1) for w < —85. The critical values of w characterizing the border of subluminal and 
supcrluminal regimes decrease for smaller n. In Fig. [5] the border line that marginally reaches the value c^ = 1 is 
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Figure 4: Evolution of the propagation speed squared c; 
Lo — —10, and (c) n = 3, uj = —1. The present epoch {z 



versus the redshift z for the cases (a) n — 2, lj = ■ 
- 0) is identified as Qde = 0.72, flm = 0.28, and Qr 



-188. 



(b) n = 2.5, 
< 10"^ 



plotted as a dotted curve. When n = 2 the avoidance of the superluminal propagation corresponds to uj > —188, see 
the case (a) in Fig. 01 Note that this condition is consistent with the result found in Ref. [26| . 

If n is close to 3, the temporal superluminal propagation can be avoided only for the values of to close to 0. In 
particular, when n = 3, we find that the superluminal propagation is inevitable for a; < 0. When n = 3, Eq. (|28p 
shows that the variable y decreases as y ex e~ = a~^ unlike the cases with n 7^ 3. In the regime x <^ y, the variable x 
increases as a; oc e'^^" = a'^ " during the radiation era. Hence, even ii x ^ y initially, there is an epoch at which x and 
y become the same order. After the system reaches this epoch, the propagation speed (|56p is no longer valid. In the 
case (c) of Fig.|3]we find that c^ starts to evolve from the value 1/2 as expected, but it enters the superluminal regime 
around the end of the radiation domination. This peculiar evolution of c^ is associated with the specific evolution of 
the variables y and x for n = 3. 



CONCLUSIONS 



We have constructed modified gravitational models of dark energy starting from the general action ([T]) without a field 
potential. The presence of a non-linear self interaction ^((/))n(/)(9^(/)9^(/)) allows us to recover the General Relativistic 
behavior in the regions of high density. The functions F{(f>), B{<p)^ and ^((/i) are restricted in the power-law forms 
given in Eq. Q from the demand of obtaining dS solutions at late times. The theory with n = 2 corresponds to 
Brans-Dicke theory with the self-interaction S,{4>) oc 4>~'^ , which was recently studied in literature [2g. If ^((/)) = 0, the 
theories with the functions F{(f)) — Afpj(0/Afpi)'^~" and B{(f>) — a;(0/Mpi)^~" are equivalent to Brans-Dicke theory 
by introducing a new field x = F((j)). However, except for n = 2, the presence of the field self interaction term does 
not allow us to express theories with the functions (|9]) as Brans Dicke theory plus the term £,{x)Ox{d'^x9tiX) by such 
a field redefinition. 

We have derived the dynamical equations ([T5)) - (pn)) to discuss the background cosmological evolution for the theories 
with the functions ^ . The evolution of the dimensionless variables x and y during radiation and matter eras can be 
analytically estimated, which is consistent with the results obtained by numerical integrations. The presence of these 
epochs demands the condition n < 3. The variables x and y at the dS point, denoted as Xds and j/ds respectively, 
are determined by solving Eqs. pUj) and (|5^ for given n and uj. We also studied the stability of the dS point by 
considering the evolution of homogenous curvature perturbations and found that the stability condition is given by 
(n - 3)a;ds < 3. 

The viable model parameter space can be restricted further by studying conditions for the avoidance of ghosts and 
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instabilities. The no-ghost condition corresponds to Qs > 0, where Qs is defined in Eq. (|47|) . Provided F(c/)) > 0, we 
require that a; > to avoid ghosts during the cosmological evolution from the radiation era to the dS epoch. We note 
that the case n = 2 is special because of the divergent behavior of the quantity Qg at the dS point. Interestingly the 
stability of the dS point is automatically satisfied for n < 3 and Xds > 0. 

The instability of perturbations can be avoided for c^ > 0, where the propagation speed squared Cg is defined in 
Eq. ([50| . At the dS point c^ reduces to Eq. (|36l). which is positive for n > 2 under the conditions n < 3 and Xds > 0. 
Hence the viable parameter region of n is constrained to be 2 < n < 3. For the existence of the dS point the parameter 
u! is restricted in the range lo < —n{n — 3)^ from Eq. (jlOp . When n ^ 3 the evolution of c^ during radiation and 
matter dominated epochs can estimated as Eqs. (j54p and (|55p respectively, which remain subluminal. During the 
transition from the matter era to the dS epoch the propagation speed can be superluminal, depending on the values 
of w and n. The avoidance of the temporal superluminal propagation gives the lower bound on cj for each value of n. 
In Fig.[2]we have plotted the viable model parameter space and the region of the subluminal propagation in the (n, uj) 
plane. When n = 3 it is difficult to avoid the appearance of the superluminal mode during the matter era because of 
the specific evolution of the variables x and y. 

It will be of interest to study the evolution of matter density perturbations and gravitational potentials to confront 
our theories with the observations of large scale structure, cosmic microwave background, and weak lensing. We leave 
the detailed analysis of cosmological perturbations for future work. 
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Appendix A: General de Sitter solutions 

We will look for the existence of general dS solutions. Let us assume the following form for (/> on the dS background 
(H = Has)- 

^ = HA^cpx{(t>). (Al) 

Once a; is a given function of 4>, Eq. (|A1[) can be easily solved as 

d0 



/ 

J d)n 



Hds{t-to). (A2) 



^ ^ :^z-2 ' (A4) 



where 4>q and H^a, are constants. Taking the time derivative of Eq. (|A1[) . we find 

(f) ^ Hl^cj) X {x + (f) X ^ff,) . (A3) 

We use the Friedmann equation ([6|) to find B, that is 

2(3i^ - 6^0^ x^ +_3£;0 (l)x + l^ (j)'^ x'^) 
(/)2 x"^ 

where ^ = H^gC Combing this with Eq. ^, it follows that 

- _ 6F + 5F <l)x + f 00 0^ x^ + f (j)x^ + f 0^ X x,0 

^" 2<f>^x^{x + <f>x,^ + 3) ■ ^ ' 

After giving x and F as functions of (j), one finds ^ and B from Eqs. (jA4p and (jASp . Thus we have proved the 
existence of dS solutions in terms of two free functions, x{(j)) and F((f>). Since the field is known as a function of t 
from Eq. (jA2p . the evolution of F, B, ^, Qs, and c^ at the dS point is also determined. In particular, the evolution 
of the time-dependent homogeneous curvature perturbations about the dS background is given by 



TZ — CI + C2 d(j) 



exp{-3/^d0[^x(0)]-i} 
QS)^xi4>) 



(A6) 



If the second term on the r.h.s. of Eq. (jA6|) decays and/or oscillates with time, then the dS solution can be regarded 
as stable. 
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Appendix B: Theories with A < 



If A < 0, then the variable y is negative from the definition ([15)) . For the existence of the matter fixed point ([25)) 
we require that n > 3. The no-ghost condition ()49p can be satisfied for a; < 0. From Eq. ([TSf we have A < for 
—3 < Xds < and n > 3. Meanwhile, when n > 3, the propagation speed squared ([5S)) is positive only in the region 
Xds < —2(4 + v^), which is not compatible with the condition —3 < Xds < 0. 

Hence the theories with A < are not viable. 
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